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Abstract 

Electron bunches with charge densities p of the order of 10 2 to 10 3 [nC/cm 3 ], energies 
between 20. and 100. [MeV], peak current > 100 [A], bunch lengths between 0.3 
and 1.8 [cm], and bunch charge of 2.0 to 20. [nC] are relevant to the design of Free 
Electron Lasers and future linear colliders. In this paper we present the results of 
numerical simulations performed with a particle in a cell (pic) code of an electron 
bunch in a drift tube. The electron bunch has cylindrical symmetry with the z-axis 
oriented in the direction of motion. The charge density distribution is constant in 
the radial and Gaussian in the longitudinal direction, respectively. The electron 
bunch experiences both a radial pinch in the middle of the pulse, corresponding to 
the peak electron density, and a significant growth of the correlated emittance. This 
behavior is explained, and an approximate scaling law is identified. Comparisons of 
the results from the pic and PARMELA codes are presented. 



1 Introduction 

New designs of Free Electron Lasers and high-brightness colliders require the 
modeling of electron bunches with higher charge densities than have previously 
been studied. In this high-charge regime, self-fields must be handled consis- 
tently. "Particle-in-a-cell" , or pic, codes must therefore be employed instead 
of the commonly-used "particle-pusher" codes, such as PARMELA, which do 
not account completely for self-fields. We studied the simple case of a bunch 
traveling along a drift tube in order to gain a better understanding of the 
self-field effects in this high-charge regime. We have completed a detailed se- 
ries of calculations using a pic code (detailed below) on a range of energy 
and electron charge where significant disagreement between pic and "particle- 
pusher" codes has been found. Our results display rather severe self-pinching 
of the bunch under certain circumstances, implying the existence of a radially 
inward-directed force which we explain analytically. 
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The outline of the paper is as follows: Section 2 contains a discussion of the 
forces acting on an electron bunch, the origin of the radial pinching, and an 
estimate of the initial space charge potential. Section 3 contains numerical 
results of the pic code, a discussion of the scaling parameters, and comparison 
between the pic and PARMELA codes results. Section 4 contains a summary 
of the paper and conclusions. 



2 Radial pinching force and space charge potential 

It is well known [1,2] that for an infinitely long relativistic electron beam, 
traveling in the positive z-direction and having constant electron density along 
z, the outward-directed radial space charge force is nearly balanced by the 
inward-directed Lorentz force. The electron beam is thus subjected to a small 
net defocusing effect. In this section we show, for the case of a finite length 
electron bunch, the origin of a radial force which is directed inward and which 
generates a pinch in the electron bunch envelope as seen in Figure 1. The 
motion of the particles is in the positive z-direction, cylindrical symmetry is 
assumed, and the relative longitudinal particle position in the bunch is denoted 
by (, where ( = z — z c , and z c is the z-position of the center of the bunch. 



(a) (b) 




Fig. 1. Electron bunch and bunch envelope at drift positions (a) z c = [cm] (drift 
entrance) and (b) z c = 30. [cm] (drift exit). The pulse envelope is projected on the 
bottom panel. Bunch parameters: Q = 2.0 [nC], a = 0.1 [cm], a = 0.1274 [cm], 
b = 1.5 [cm], E = 20.0 [MeV], 70 = 40.138 . 

We assume, following [3], that all the electrons travel along z with uniform 
velocity v = fioc, no transverse velocity component is present, and the density 
p(r, z) is given by 

Q/(vra 2 ) f(s) = nj(s) : r < a 
: r > a, 

where s = z — /3 ct, f(s) = e~( s /°") 2 / 2 /[(27r) 1 / 2 (j] is the normalized longitudinal 
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p(r,s(z,t)) = < 



density distribution, Q is the pulse charge with appropriate sign, and a is the 
electron bunch radius. We further assume that the z- and t-dependencies of 
the electric field E and magnetic field B in Maxwell's equations are combined 
in the variable s. This ansatz is appropriate when the electrons have a common 
longitudinal velocity (3q and no transverse motion. This assumption is violated, 
as shown by the numerical simulations, after the bunch drifts a significant 
distance. It is, however, useful to establish the conditions for the onset of the 
radial pinching force. The region in which we will solve Maxwell's equations 
is the inside of a perfectly conductive drift tube of radius b. 

In view of the symmetry of the problem, only the E®, E®, and B$ components 
of the fields are present. The superscript ()° is used to indicate that the fields 
are computed in the absence of any transverse velocity. 

The Fourier transform of the electric field is 

£(k, r) = (27T)- 1 J E(r, s) e~ iks ds, (2) 



and similar definitions hold for the Fourier transforms B(k, r) and T{k) of 
B(r, s) and f(s), respectively. For our selection of /(s), a Gaussian distribution 
with standard deviation a is assumed, and J-(k) = e~( feCT ) I 2 . 

Upon substitution into Maxwell's equations, £® and Bg can be expressed as 
functions of £® only, and inside the perfectly conductive drift tube the results 
are as follows: 



_i/c<9£^ 
_ikbd£l 

where q = — i/c/70 and 7q = (1 — /3g) _1 . In the region inside the electron bunch, 
where r < a, the differential equation satisfied by £® is 



and in the region outside the electron bunch, where a < r < b, 

+ q 2 £° z = 0. (6) 




The solution of Eqs. (5) and (6) must be finite, continuos everywhere, vanish 
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at the conductive wall of the drift tube (r = b) and insure the continuity of 
£® at the edge of the pulse (r = a). The solution is 



£°,(k,r) 



CMqr) + 47m T{k)/(ik) 
[d - 27r 2 n gar 1 (ga) ^(fc)/(iA;)] J (gr) 
+[27r 2 n ?aJi( g a) F(k)/(ik)] F (gr) 



r < a 



r > a, 



(7) 



where J n and F„ are the Bessel functions of first and second kind [4]. The 
integration constant C\ is 



Ci = 2n 2 n qa x{o) ^(k)/(ik), where 
X (q) = \Vi(qa) J (qb) - J^qa) Y (qb)} / J (qb). 



(8) 
(9) 



2. 1 Calculation of the radial force 



The zeroth order calculation of the radial force, _F S °, neglects the electron radial 
motion induced by the space charge and Lorentz force, i.e. the electron bunch 
moves as a solid body. F s ° is the sum of the radial space charge force and the 
Lorentz force arising from the poloidal magnetic field B§ and the longitudinal 
velocity of the electrons: 

F° = e(E° r -(3B° e ), (10) 



where e is the electronic charge with appropriate sign. With the help of Eqs. 
(3), (4), and (7) (for r < a), we find for the radial force 



F s ° = -(l-/5 2 )^ f dke^-^^xiikh^iikrho) 

— oo 

= (1 - (11) 

It is simple to verify that the integral is always real and negative, thus the 
radial force is directed outward and corresponds to a defocusing of the electron 
pulse. This result is a simple generalization of the one obtained in [1,2] for an 
infinitely long electron beam. 

The first order correction to the radial force includes the effect of the electron 
radial motion. Inspection of the transverse phase space plots, (x',x), where 
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x> — fix/ fiz, obtained from the numerical calculations reveals that the electrons 
acquire a significant radial motion (see Figure 6). Let 

Be = B° d + Bl (12) 



where B\ is the correction term arising from a non- vanishing radial velocity 
of the electrons. B\ is obtained from Ampere's law 



whose projection along the radial direction is 

+ H Jr . (14) 

az c at c 



We substitute Eq. (12) into (14) and observe that the homogenous part of Eq. 
(14) is satisfied by Bg, and the result is 

= -^J, (15) 

az c 



The radial momentum equation for the electrons is -^ t (m^v r ) = eE r (l — (3 Q ). 
The numerical calculations indicate that the electron radial excursion and its 
energy change is quite limited in the initial phases of the pinch. Hence, the 
radial momentum equation can be integrated at the onset of the constriction 
with the assumption that 70 and E r are constant along the electron trajectory. 
The resulting j r is 

3r = ~ 6 1 ~ f3 ° p(z,r)E%z, (16) 
7 mc P 



where the relationship = — (3qc-^ has been used. We substitute Eq. (16) 
into (15) and integrate, and the result, at the center of the electron bunch 
(C = 0), is 

^ = (J)>, (17) 



where, at z c = 0, the electron bunch is at the beginning of the drift tube, and 



'2 nI A a 2 a 7/? \ 1/2 
Wc 1-/3 J 
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9 \ V 2 

^V/J , (18) 



V 2g c 

where I a = mc 3 /e is the Alfven current with value 1.7 x 10 4 [Amps], or 
3.12 x 10 13 [StatAmps]. 

We substitute Eqs. (18) and (12) into the Lorentz force and obtain the first 
order correction to Eq. (11), the electron radial force: 



(19) 



Eq. (19) indicates that, at sufficiently large z c values, the radial force will be 
turned inward. The expression for z s shows that high charge and low energy 
beams will become more constricted than their counterparts with lower charge 
and higher energies. For an electron bunch with Gaussian longitudinal charge 
density distribution, such as the one considered in this paper, E® peaks at the 
center of the pulse. Consequently, the maximum constriction will also occur 
at the center of the pulse as has been observed in the numerical calculations. 



2.2 Calculation of the space charge potential 



An estimate of the initial space charge potential at the beginning of the drift 
region has been made via the following procedure. The space charge potential 
$' is first computed in the electron rest frame K' which moves along z with 
velocity v = fioc relative to the laboratory frame K. In K', the electrons are 
stationary, the electron relative motion is neglected, the electric potential is a 
well-defined quantity, and it satisfies the Poisson equation. $' is then trans- 
formed back into the laboratory frame K, where, from energy conservation, it 
is translated to an equivalent A7. The details of the procedure are as follows: 

C C 
$'(C',r) = - J d V E' z ( V ,r) = <5>'((,r) = - 1 J d V E z ( V ,r), (20) 

—00 —00 



where &((,r) transforms as the fourth component of a vector in the four- 
dimensional space-time space. Thus, we obtain 

C 

$(C,r) = - 7 2 J d V E z ( V ,r). (21) 

—00 
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Energy conservation gives 



7 mc 




— oo 



E z (r),r) is obtained from Eq. (7) (for r < a). Thus, 



mc z a J v > 

— oo 

x [7rx(?)^o(?0 + 2 7 /(ia£;)] , 



(23) 



where Co is arbitrary, provided it satisfies the condition ( » a. Finally, we 
compute the normalized momentum z-component 



7 & = ( 7 o + A 7 (C,r)) 2 -l 



1/2 



(24) 



In the derivation of Eq. (24), the electrons are assumed to be strictly stationary 
in K'. When significant radial motion is present, the results of this analysis 
will be inaccurate as is observed in the numerical results. 



3 Numerical simulation results 



The numerical calculations of a cylindrical, axisymmetric electron bunch in 
a drift tube with perfect wall conductivity were performed with the "particle 
in a cell" code ARGUS [5]. ARGUS is a fully three-dimensional and time- 
dependent (x, y, z, t) solver for Maxwell's equations. We used the "pic" solver 
option. 

We also compared ARGUS results to the calculations from the standard 
code PARMELA [6]. PARMELA was originally developed at Los Alamos Na- 
tional Laboratory and is widely used within the accelerator community. In 
PARMELA, the forces between the electrons are computed in the electron 
rest frame where all the electron relative motion is neglected. 



3.1 ARGUS results 



The electron bunch is generated at z — 0, and the electron density conforms to 
Eq. (1). The total electron bunch length in our simulation is 5.0 cr's. a varies 
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from 0.1274 to 0.7644 [cm]. The bunch radius a varies from 0.1 to 0.3 [cm], 
the electron energy E from 20. to 100. [MeV], and the charge Q from 2.0 to 
18.0 [nC]. The electrons are injected with uniform f3 z , and (5 x — (5 y — 0. 

In the code the electrons are represented by "macro particles". The density 
of the macro particles is uniform, but the charge is weighted to reproduce 
the appropriate charge density distribution. The drift tube, as represented in 
ARGUS, is shown in Figure 2. 




Fig. 2. Drift tube geometry in ARGUS. The positive z-axis is oriented in the direc- 
tion of the bunch motion. 

The drift tube radius is 1.5 [cm] and the length is 30.0 [cm]. The number of 
grid points along x and y is 57, and the number along z varies from 192 to 
571. The spacing in x and y is not uniform: the smaller spacing (0.02 [cm]) is 
used in proximity of the z-axis of the tube, and the larger spacing (0.163 [cm]) 
in the vicinity of the wall. The spacing along z is constant. The time step, 
0.35 [ps], is selected such that the Courant condition is satisfied everywhere. 
This condition is v p At < Ax min /2, where v p is the particle velocity, At is the 
time step, and Ax min is the smallest grid spacing. 

In Figure 3 the evolution of the (x, () bunch envelope at different positions 
in the drift is shown. In view of the cylindrical symmetry of the electron 
bunch, the (x, () envelope is identical to the (y, () one. At z c = 0. the cross 
section is uniform, in the interval 10. < z c < 20. [cm] a small pinch develops, 
and at z c = 30. [cm] a fully evolved pinch has been established. The largest 
pinch occurs in the middle of the bunch in correspondence to the maximum 
radial force as also seen in Eq. (19). A close inspection of the bunch envelope 
reveals a tapering of the rear half. The wake fields from the front half may be 
responsible for this behavior. 
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Fig. 3. Evolution of the bunch envelope. Bunch parameters: Q = 2.0 [nC], a = 0.1 
[cm], a = 0.1274 [cm], b = 1.5 [cm], E = 20. [MeV], 70 = 40.138 . 

In Figure 4 the evolution of the longitudinal phase space, 7/5 2 vs (, is pre- 
sented for the same case as in Figure 3. In the inset of Figure 4(a) we have 
overlapped, with an expanded scale, the results of Eq. (24) with the ARGUS 
calculations for purpose of comparison. Here, the lower curve corresponds to 
r = and the upper one to r = a. The agreement is quite satisfactory in 
view of the simplifying assumptions made in our analytical formulation. As 
the pulse drifts in the tube, the relative electron motion increases, and a more 
complex longitudinal phase space evolves. 
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Fig. 4. Evolution of the longitudinal phase space at drift positions (a) z c = 0, (b) 
z c = 10. [cm], (c) z c = 20. [cm], and (d) z c = 30. [cm]. Bunch parameters: Q = 2.0 
[nC], a = 0.1 [cm], a = 0.1274 [cm], b = 1.5 [cm], E = 20. [MeV], 7o = 40.138 . 
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In Figure 5 the phase space (x',x,z) is shown at z c = and z c = 30. [cm]. 
At z c = the phase space is planar as expected from the initial conditions 
inposed on the momemtum of the electrons. At z c = 30. [cm] the phase space 
has become twisted. A projection of the intersection of the phase space with 
the C = plane, where the largest pinch occurs, is displayed on the rear panel 
of the figure. A significant inward motion of the electrons is present. 




Fig. 5. Evolution of the (x', x, Q phase space at drift positions (a) z c = and 
(b) z c = 30. [cm]. The pulse envelope is projected on the bottom panel. Bunch 
parameters: Q = 2.0 [nC], a = 0.1 [cm], a = 0.1274 [cm], b = 1.5 [cm], E = 20. 
[Me V], 70 = 40.138 . 

In Figure 6 the intersections of the phase space (x', x, z) with planes at con- 
stant ( are shown when the pulse is at the end of the drift tube, z c = 30. 
[cm]. Figure 6(a) and 6(b) correspond to the tail and the middle of the bunch, 
respectively. The entire (x',x) phase space is shown in Figure 12(a). 
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Fig. 6. Intersections of the (x',x,C) phase space at drift position z c = 30. [cm]. 
Intersections are at (a) z = —0.315 [cm], (b) z = —0.030 [cm], (c) z = +0.094 [cm], 
and (d) z = +0.315 [cm]. Bunch parameters: Q = 2.0 [nC], a = 0.1 [cm], a = 0.1274 
[cm], b = 1.5 [cm], E = 20. [MeV], 7o = 40.138 . 
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It is convenient to introduce the concept of the fractional rms "slice" emit- 
tance. The total normalized rms emittance of the pulse is defined as 



-AT 



= 7/9 



7T 



X 



x' 2 



1/2 



(25) 



where the averages indicated by the angle brackets, (), extend over all the 
pulse particles. The fractional slice emittance is defined as 



£ rms,f(C, AC) = 7/3 



7T 



x 2 ) 



>2 



X ' 



AC 



(xx 1 ) 



2 

AC 



1/2' 



N 

erms ' 



(26) 



where the averages indicated by the angle brackets, () A ^, are limited to parti- 
cles within the slice, i.e. with ^-coordinate between ( and ( + A(. 

In Figure 7(a) the pulse current profile at z c = and z c = 30. [cm] is shown. 
The current profiles are nearly identical, as expected, because of the limited 
relative longitudinal motion among the particles. In Figure 7(b) the fractional 
slice emittance as a function of ( for different position z c along the drift tube 
is shown. Observe that the fractional slice emittance is always significantly 
smaller than the total pulse emittance of 18.5 [ir mm-mrad] at the end of the 
drift tube. This indicates that the pulse emittance is strongly correlated. 
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Fig. 7. (a) Bunch current and (b) fractional slice emittance vs position in the bunch 
at different locations in the drift tube. All fractional slice emittances in (b) are 
normalized to 18.5 [tt mm-mrad]. Legend: dark solid line, z c = 0.66 [cm]; dark 
broken line, z c = 10. [cm]; light broken line, z c = 20. [cm]; light solid line, zc = 30. 
[cm]. Bunch parameters: Q = 2.0 [nC], a = 0.1 [cm], a = 0.1274 [cm], b = 1.5 [cm], 
E = 20. [MeV], 70 = 40.138 . Curve fluctuations are due to discretization noise. 

In Figure 8 typical longitudinal (C,^) bunch envelopes at zc = 30. [cm] are 
presented, corresponding to varying energies, bunch lengths, and charges. The 
pinch decreases at higher energies and increases with bunch length for constant 
charge density. 
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Fig. 8. Various longitudinal (C,%) bunch envelopes at drift location z c = 30. [cm]. 
For (a)-(d), a = 0.1 [cm]. Bunch parameters: 

(a) Q = 2.0 [nC], a = 0.1274 [cm], b = 0.25 [cm], E = 50. [MeV], 7o = 98.84 

(b) Q = 2.0 [nC], a = 0.1274 [cm], b = 0.25 [cm], E = 100. [MeV], 70 = 196.69 

(c) Q = 2.0 [nC], a = 0.7639 [cm], b = 0.25 [cm], E = 20. [MeV], 70 = 40.13 

(d) Q = 6.0 [nC], a = 0.3822 [cm], b = 1.5 [cm], E = 20. [MeV], 7o = 40.13. 

The evolution of the pinch ratio as a function of the position in the drift tube is 
shown in Figure 9(a). Here, r min /a, the ratio of the minimum bunch envelope 
radius divided by the initial bunch radius, is plotted versus the position in the 
drift tube for three different values of the bunch energy. Higher energy bunches 
undergo less pinching than lower energy ones. In Figure 9(b), r m i n /a is plotted 
versus the normalized bunch center position, z c /z s . The three separate curves 
of Figure 9(a) nearly collapse into a single one, thus showing that Eq. (18) 
correctly captures the energy scaling of the pinch ratio. 
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Fig. 9. Pinch ratio vs (a) drift position z c and (b) normalized drift position z c j z s 
in the drift tube for various initial energies. Common bunch parameters: Q = 2.0 
[nC], a = 0.1 [cm], a = 0.1274 [cm], b = 0.25 [cm]. 
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In Figure 10 the evolution of the pinch ratio is shown as a function of z c j z s 
for pulses of different diameters, same lengths and different charges selected to 
maintain a constant electron density. As in the previous case, the pinch ratio 
is well parameterized by z c jz s . 
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Fig. 10. Pinch ratio vs normalized drift position zc/zs in the drift tube for various 
charges and initial bunch radii. Common bunch parameters: a = 0.1274 [cm], b = 1.5 
[cm], E = 20. [MeV], 70 = 40.138. 

Finally, in Figure 11 the evolution of the pinch ratio is shown as a function 
of z c /z s for bunches of different lengths, identical transverse cross sections 
and different charges selected to maintain a constant electron density. Unlike 
the previous cases, the pinch ratio is not well parameterized by z c /z s . The 
probable explanation for the inability to parameterize the present case is that 
the wake fields of the particles in the leading edge of the bunch are not properly 
accounted for in our simple theory. 
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Fig. 11. Pinch ratio vs position z c j z s in the drift tube for various charges and bunch 
lengths. Common bunch parameters: a = 0.1 [cm], b = 1.5 [cm], E = 20. [MeV], 
70 = 40.138. 
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3.2 Comparison of ARGUS and PARMELA results 



ARGUS and PARMELA yield quite different results for identical bunch con- 
ditions. In Figure 12 the transverse phase space, as computed by ARGUS and 
PARMELA at z c = 30. [cm], is shown. The PARMELA phase space reflects 
only a mild defocusing effect due to the space charge. The ARGUS counterpart 
combines the focusing and defocusing due to the pinching force. 
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Fig. 12. (a) ARGUS and (b) PARMELA results of transverse (x, x') phase space 
at drift position z c = 30. [cm]. Bunch parameters: Q = 2.0 [nC], a = 0.1 [cm], 
a = 0.1274 [cm], b= 1.5 [cm], E = 20. [MeV], 70 = 40.138. 

In Figure 13 the longitudinal phase space, as computed by ARGUS and 
PARMELA at z c = 30. [cm], is shown. Here, PARMELA cannot reproduce 
the complexity of the longitudinal momentum distribution within the bunch. 
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Fig. 13. (a) ARGUS and (b) PARMELA results of longitudinal {<yf3 z ,() phase space 
at drift position z c = 30. [cm]. Bunch parameters: Q = 2.0 [nC], a = 0.1 [cm], 
a = 0.1274 [cm], b = 1.5 [cm], E = 20. [MeV], 7o = 40.138. 



14 



Finally, in Figure 14 the longitudinal ((,x) envelopes at z c = 30. [cm], as 
computed by ARGUS and PARMELA, are presented. No pinching of the 
bunch is exhibited by the PARMELA calculations. 
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Fig. 14. (a) ARGUS and (b) PARMELA results of longitudinal (C, x) bunch en- 
velopes at drfit position z c = 30. [cm]. Bunch parameters: Q = 2.0 [nC], a = 0.1 
[cm], a = 0.1274 [cm], b = 1.5 [cm], E = 20. [MeV], 7o = 40.138. 

4 Summary and Conclusions 

We have presented the results of extensive numerical simulations of electron 
bunches with charge densities p of the order of 10 2 to 10 3 [nC/cm 3 ] and energies 
between 20. and 100. [MeV]. The results indicate the presence of a strong 
pinch in the middle of an electron bunch with a Gaussian longitudinal electron 
density distribution. The pinching force scales approximately as 7~ 3 / 2 and r 1 / 2 . 
This force generates an increase in the correlated bunch emittance, and the 
space charge depression is also affected. A simplified analysis to explain these 
results has been described. The results are in significant disagreement with 
the ones obtained from the standard accelerator code PARMELA. This is 
attributed to the neglect of the electron relative motion in the space charge 
calculations performed in PARMELA. 

The pinch effects we have uncovered will play a significant role in the pho- 
toinjectors used in Free Electron Laser and other high-brightness accelerators. 
Here, the electron pulse has very low energy, and this will limit the allowed 
charge density in order to avoid undue pinching and space charge effects. The 
break up of the electron bunch at the photoinjector has been observed [7,8] 
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and attributed solely to the longitudinal component of the space charge force. 
The pinching force discussed in this paper may also play a role. In addition, 
the design of magnetic pulse compressors will need to account for the pinch 
which may become significant in the last stage of the compression. 

Further analytical studies and numerical simulations are necessary to better 
understand (a) the evolution of the pinch as z c — > oo and (b) the effect of the 
longitudinal electron density distributions. Preliminary indications are that 
uniform longitudinal charge density distributions become pinched at multiple 
positions in the bunch. The pinch period depends upon the pulse parame- 
ters. More importantly, experimental observation of the pinch is needed to 
completely validate these results. 
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